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We study the thermal and non-thermal steady state scaling functions and the steady-state dy¬ 
namics of a model of local quantum criticality. The model we consider, i.e. the pseudogap Kondo 
model, allows us to study the concept of effective temperatures near fully interacting as well as weak- 
coupling fixed points. In the vicinity of each fixed point we establish the existence of an effective 
temperature -different at each fixed point- such that the equilibrium fluctuation-dissipation theo¬ 
rem is recovered. Most notably, steady-state scaling functions in terms of the effective temperatures 
coincide with the equilibrium scaling functions. This result extends to higher correlation functions 
as is explicitly demonstrated for the Kondo singlet strength. The non-linear charge transport is also 
studied and analyzed in terms of the effective temperature. 
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The interest in understanding the dynamics of strongly 
correlated systems beyond the linear response regime has 
in recent years grown tremendously. 

The quantum dynamics in adiabatically isolated opti¬ 
cal traps has been successfully modeled using powerful 
numerical schemes El- In open systems mainly di¬ 
agrammatic techniques on the Schwinger-Keldysh con¬ 
tour have been employed. For nanostructured systems 
several techniques exist to describe the ensuing out-of¬ 
equilibrium properties. These approaches, however, are 
either perturbative in nature [3], centered around high 
temperatures and short times SHU, or approximate the 
continuous baths by discrete Wilson chains [SHU. The 
situation might be simpler for non-linear dynamics that 
arises in the vicinity of a quantum critical point (QCP), 
where a vanishing energy scale leads to scaling and uni¬ 
versality [T2UT9] . 

For the dynamics near classical continuous phase transi¬ 
tions a well-established theoretical framework exists, ty¬ 
ing the dynamics to the statics and the conserved quan¬ 
tities m- In addition, the concept of effective temper¬ 
ature (T efI ) was established as an useful notion for the 
relaxational dynamics of classical critical systems E3- 
123] , although it appears somewhat less useful for fully 
interacting critical points [23] . T off is commonly defined 
by extending the equilibrium fluctuation-dissipation the¬ 
orem to the non-linear regime. The existence of effective 
temperatures in quantum systems was recently investi¬ 
gated |TgI[25H2g| . For a recent review see [25] ■ Ti compar¬ 
ison to classical criticality, at a QCP, dynamics already 


enters at the equilibrium level. For a QCP that can be 
described by a Ginzburg-Landau-Wilson functional in el¬ 
evated dimensions, it was found that the voltage-driven 
transition is in the universality class of the associated 
thermal classical model with voltage acting as T eff m- 
Unconventional QCPs in contrast are not described solely 
in terms of an order parameter functional |291150] . 

In this letter we address the following general ques¬ 
tions within a model system of unconventional quantum 
criticality: Is the existence of T efr tied to dynamical (or 
w/T-)scaling? Does T afr have meaning for higher correla¬ 
tion functions? How unique is T efr at a given fixed point 
once boundary conditions have been specified? Can crit¬ 
ical scaling functions be expressed through T off and if so, 
how do these scaling functions relate to the equilibrium 
scaling functions? The model system is the pseudogap 
Kondo (pKM) model that describes a quantum spin anti- 
ferromagnetically coupled to a conduction electron bath 
possessing a pseudogap near its Fermi energy, character¬ 
ized by a powerlaw exponent. Depending on the coupling 
strength, the quantum spin is either screened or remains 
free in the zero temperature (T) limit. The two phases 
are separated by a critical point dispaying critical Kondo 
destruction, see FigjTj The pKM has been invoked to 
describe non-magnetic impurities in the cuprate super¬ 
conductors m and point-defects in graphene [32!. It 
underlies the pseudogap free moment phases occurring 
in certain disordered metals [33} and can also be realized 
in double quantum-dot systems [34]. The quantum crit¬ 
ical properties of the pKM in equilibrium have been ad- 
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Figure 1. (a) Sketch of the model: a spin interacts with two 
fermionic leads which are characterized by their respective 
density of states p~ L / R (oj) and chemical potential Rl/r- (b) 
Phase diagram of the multichannel pKM with gap exponent 
r < r max : A QCP (C) separates the multichannel Kondo fixed 
point (MCK) from the (weak-coupling) local moment fixed 
point (LM). 


dressed in [4511451 . Our main findings are that the steady 
state dynamic spin susceptibility, the conductance, and 
the Kondo-singlet strength, a 4-point correlator, repro¬ 
duce their equilibrium behavior in the scaling regimes of 
the fixed points of the model when expressed in terms of 
a fixed-point specific T eft . 

The model. We consider a pKM with a density of 
states that vanishes in a power-law fashion with ex¬ 
ponent 0 < r < 1 at their respective Fermi level, 
p~l(cj) ~ |cj| r 0(D — |w|), with half-bandwidth D. Here, 
l = L,R labels the two leads, see Fig{l])a). In the 
multichannel version of the model the spin degree of 
freedom ( S ) is generalized from SU( 2) to SU(N) and 
the fermionic excitations (c) of the leads transform un¬ 
der the fundamental representation of SU(N ) x SU(M) 
with N spin and M charge channels. At T = 0 and 
r < r max < 1, a critical point (C) separates a multichan¬ 
nel Kondo (MCK)-screened phase from a local moment 
(LM) phase at a critical value J c of the exchange cou¬ 
pling J > 0, see Figfljb). The characterization of the 
phases and the leading power law exponents of observ¬ 
ables of the pKM have been obtained by perturbative 
RG, large -N methods, and NRG [55H581PUJI145] . Within 
the large-IV approach, at T = 0, scaling arguments are 
able to predict the critical exponents of dynamical ob¬ 
servables [391 [46] . Non-equilibrium steady-states (NESS) 
are obtained by applying a time-independent bias voltage 
V = (pl — Pr)/ |e|, where /q is the chemical potential 
of lead l, see Fig. [TJa) . As T characterizes the fermionic 
reservoirs, it remains well-defined even for V ^ 0. 

A similar setup has been considered in a perturbative 
RG-like study adapted to the NESS condition [47 . This 
model has also been invoked in a variational study of the 
dynamics following a local quench where it was found 
that quenches in the Kondo phase thermalize while this 
in not the case for quenches across the QCP into the LM 
regime |48|. The system is described by the Hamiltonian 

H = S £ P lC paal c P<x<Tl + Jll'S.S a .U', (1) 

p,aal IV ot 




Figure 2. (a) x' (0) _1 vs J for different T. (b) <j> B vs J for 
different T. The T — 0 curve is approached from below in 
the MCK and from above in the LM phases, (c) Scaling 
T/T b s vs V/T at fixed points LM, C, and MCK: T bS ~ V 
for V T. (d) FDR7 1 (w) vs w/T e fj near fixed point C, 
shown for V/D = 10 -2 ,10 —3 ,10 -4 ,10 -5 ,10 -6 . The grey line 
is FDR~ 1 (cu) = tanh(/3w/2). 


where a = 1 and a = 1 are, respec¬ 

tively, the SU (A^-spin and SU(M)- channel indices, l 
labels the leads and p is a momentum index. The co¬ 
tunneling term (49] in Eq. ([!]) contains the local opera¬ 
tors <•«' = ^Epp'aa' with * the fun¬ 

damental representation of SU(N) and n c is the num¬ 
ber of fermionic single-particle states. In a totally anti¬ 
symmetric representation, one can decompose the spin 
operator as S a(7 > = — qS aa ', where q is subject to 

the constraint Q = ftf o- = 9 N and the /t, f a , obey 

fermionic commutation relations. 

We employ a dynamical large-N limit [59} [50] , suitably 
generalized to the Keldysh contour mm while keeping 
q = % and k = M/N constant. This results in 

S>* < (t) = *G>- < (t)G<->(-t) (2) 

E >’ < (t) = -iKG%’ < (t)G>’ < (t) (3) 

-iGf (0) = q, (4) 

where G/ is the pseudofermion propagator and Gb is 
the propagator of a bosonic Hubbard-Stratonovich de¬ 
coupling field. E f (Eb) is the proper selfenergy of G/ 
(Gb) and is related to it via the Dyson equation [ 5T] . 
We assume that the exchange interaction originates from 
an Anderson-type model via a Schrieffer-Wolff transfor¬ 
mation, so that a single coupling constant J = Jl + Jr 
emerges m- 

For details on the numerics see m- In equilibrium, our 
approach yields dynamical scaling functions that coincide 
with those obtained from quantum Monte-Carlo [44]. 

Observables. A possible order parameter for the transi¬ 
tion from the overscreened Kondo to local-moment phase 
is given by limr^o Tx(u = 0, T), where x(w,T) is the 
Fourier transform of the local (impurity) spin-spin cor¬ 
relation function x(t ~ t 7 ), see Fig. !a). We work on 
the Keldysh contour where the lesser and greater com¬ 
ponents are defined in the usual way as y > (t — t') = 
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Figure 3. Scaling of observables with T e g at different fixed 
points for the values of V as in Fig. [ 2 J d): (a) x! (0) _1 vs T e fr; 
(b) 9i n & In x" (w) vs w/T e fr; (c) c fi 3 vs T e g. For each fixed 
point, the equilibrium scaling form (grey curves) is compared 
with the same quantity under non-equilibrium conditions and 
T substituted by T e g. 


with t G 7<- and t' € 7-*- and 
x< (i - t') = E« (S a (O s a (t)), with t G 7 -v and 
t' G 7<_ so that x^ 0) = © (i) [X > (*) - X < (i)] and 
X" 4 = X R + X < — X > • Here, is the forward (back¬ 

ward) branch of the Keldysh contour, respectively. 

We also consider the “singlet-strength” cf> s , defined 
through the Kondo term contribution to the total 
energy of the system as ^ Ezi'X] c ( s - s c,«') = 

—Jk J Ps |52| . <f> s is a dimensionless quantity, 

which possesses a well-defined large -N limit and quan¬ 
tifies the degree of singlet formation. In terms of 
the fermionic fields, it can be written as the local- 
in-time limit of a 4-point correlator m- Its equilib¬ 
rium properties will be discussed below. The steady 
state charge current passing through each channel is 

J P = -dt(N L (t)^ /M, where N L = E pQCT c U*L c p a *L 
is the number of particles in the left lead. The out-of¬ 
equilibrium conditions considered here respect particle- 
hole symmetry which implies a vanishing energy current. 

Throughout the paper we set n = 0.3, r = 0.2, and 
q = 1/2. This results in r max = 0.412(4). Our choice of 
values for k and r ensures a finite static spin suscepti¬ 
bility x' ( w — 0) within the MCK phase as T —> 0. We 
denote the real (imaginary) part of \ r (uj) by x! (x")- 
Thermal steady-state. The equilibrium (V = 0) behav¬ 
ior of x! (w = 0 ,T) in the relaxational regime (w <C T) 
near the MCK, C, and LM fixed points is shown in 
Figj2}(a). For J < J c ~ 0.44D, i.e. in the LM phase, 
one observes Curie-like behavior at lowest temperatures 
x'(w = 0,T) oc T -1 . In the MCK phase (J > J c and 
with our choices of n and r), the T = 0 susceptibility re¬ 
mains finite. The grey lines in Figsj3j(b) show the scaling 
plots of the logarithmic derivative of x" ( w ) for different 


values of the temperature, i.e. d\ n M In x" {ui) for the dif¬ 
ferent fixed points. Note that d\ BU] In x" (w) — a x within 
the scaling region where x" ( w ) oc The values of a x 

in the quantum coherent regime (u/T 3> 1) agree with 
those obtained analytically from a T = 0 scaling ansatz 
[45 for the MCK (a x ~ 0.087) and C (a x = —0.97) 
fixed points. These results are compatible with a dy¬ 
namical scaling form x" (w, T) = T a *4> (w/T), in terms 
of an universal scaling function $ (x) possessing asymp¬ 
totic values 4> (x) oc x for x -C 1 and $ ( x ) oc x ax for 
i> 1. Thus, the scaling properties are in line with dy¬ 
namical w/T-scaling for the C and MCK fixed points. For 
the LM fixed point we find a x = — 1 and a scaling form 
x" (w) = (w/T 1+K ), indicative of a weak-coupling 

fixed point and absence of hyperscaling. These results 
will be further addressed elsewhere m- 
The singlet-strength (f> s vs. J at different T and at V = 0 
is shown in Fig[2]-(b). The numerical data at T ^ 0 sug¬ 
gest that (f) s (J, T = 0) is a continuous function of J. At 
the C fixed point we find that <f> s (J, T) as a function of 
J crosses for different values of T (for sufficiently low T). 
Non-thermal steady-states. We consider a non¬ 
equilibrium setup where the two leads, initially decou¬ 
pled from the impurity (for t < to ), are held at chemical 
potentials /i P = ~/j,r = \e\ V/2 (|e| = 1 in the follow¬ 
ing). At t = to the coupling between the leads and 
the impurity is turned on. A steady-state is reached 
by sending to —> —00 so that any transient behavior 
will already have faded away at (finite) time t. The 
NESS fluctuation-dissipation ratio (FDR) for a dynam¬ 
ical observable A(t,t') = A(t — tf) is defined through 
FDRa(w) = [A > (u) + A < {uj)\/[A > ( u) — where 

A>/< are the Fourier transforms of the greater/lesser 
components of A. At equilibrium, the fluctuation- 
dissipation theorem implies FDRx(w) = tanh(/3w/2)^ 
uniquely (with (f = ±1 for fermionic (+) and bosonic 
(-) operators). For a generic out-of-equilibrium sys¬ 
tem, the functional form of the FDR differs from the 
equilibrium one. A frequency-dependent “effective tem¬ 
perature”, l//3/g(u;), for the observable A can be de¬ 
fined by requiring that tanh [/^ (w) oj/2] ^ = FDRx(w) 
[271155] , Following Refs. [T511211125] we define T aft via 
FDR X through its asymptotic low-frequency behavior 
T~q = lim^-^o /?eff( w )- equilibrium T eff = T. On the 
other hand, a linear-in-K decoherence rate in the non¬ 
equilibrium relaxational regime near an interacting QCP 
is signaled by w/K-scaling |16) . In this case and at T = 0 
one expects T eff = cV, where c characterizes the under¬ 
lying fixed point. We thus analyze T/T eff vs V/T. Fig. 
i( c) shows the resulting T/T e g as a function of V/T for 
the different fixed points computed for different values of 
V and T. In the non-linear regime, the scaling collapse 
for T/T bB implies T off = cV, where 1/c is the amplitude of 
the scaling curve in the non-linear regime. A comparison 
of FDR” 1 with the equilibrium result for fixed point (C) 
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Figure 4. Conductance G normalised to the MCK fixed point 
conductance 2-7 tGo = 0.415. (a) G(T) vs T computed for the 
lowest non-zero value of V at different values of J (see color 
coding), (b) G vs V for fixed T. (c) G = Jv IV vs T e ff at 
different fixed points. The equilibrium form is given by the 
grey curves. 


is shown in Fig{2]-(d). Even for the LM fixed point, where 
hyperscaling is violated, T eff ~ V holds for V T, see 
Fig. §( c), top panel. It is however important to realize 
that the properties we see in terms of T efI are a property 
of the flow towards the LM fixed point. Far from equi¬ 
librium and outside any scaling regime, \ is a function of 
w, T, and V but near a fixed point %(w,T, V) develops 
a scaling form in terms of a combination of u>, T, and 
V. This then raises the question how T e g enters the scal¬ 
ing function and leads us to a remarkable result, see Fig. 
|3]-(b)-(c): The non-thermal steady-state scaling function 
of x = x(w,T,V) when scaled in terms of T eff recovers 
the equilibrium scaling function of that particular fixed 
point with T eff replacing T. This not only turns out to 
be true for \ at each of the fixed points of the model but 
also holds for </> a , a higher-order correlation function. We 
first consider the static susceptibility. Fig. [3]- (a) shows 
the equilibrium scaling forms of x! (0) _1 as a function of 
T e g for different values of T and V for the LM, C and 
MCK fixed points. The color coding reflects the values 
of T of the system. The equilibrium form (grey lines) is 
recovered even for T e g/T 1. 

A similar result can be obtained at finite w: Fig. [3} 
(b) shows the log-derivative <9i nw In x" (w) as a function 
of u)/T e g for different values of T and V for the LM, C 
and MCK fixed points. These should be compared with 
the equilibrium results, the underlying grey lines: The 
equilibrium scaling form is recovered by replacing T by 
T e ff, both for u> <C T e ff and to T e g [52]. Note that 
T e ff is defined from the FDR of x in the limit u>/T —> 0. 
Therefore, the fact that the equilibrium scaling forms 
of x! (0) and x" ( w ) are reproduced for T e g/T 1 and 
uj/T e g 1, respectively, is remarkable. Fig. [3]- (c) depicts 
tt> s as a function of T e g for different values of T and V. 
Again, the equilibrium scaling behavior (gray curves) is 
reproduced. 

Unlike x and <j> a , the conductance G depends on both 
pseudoparticle propagators Gf and Gb- One thus may 
wonder if T eff can have any meaning for G. In Figs. [4} 


(a,b) we show the conductance per channel G = Jp/V 
vs T and V respectively. In the linear response regime 
V, T -C T k (J) of the MCK phase, the current is pro¬ 
portional to the applied voltage J\> = GqV. Outside 
of the scaling regime, i.e. for F,T> Tk (J), G drops 
rapidly as V or T increase. The linear and non-linear 
current-voltage characteristics display power-law behav¬ 
ior as T, V —> 0 [THl H8]. Near C, i.e. for J = J c , the 
relation between J\> and V is still linear, (Jy> = G C U), 
however the critical conductance G c is much smaller than 
Go- Fig. |4]-(c) shows G vs T e g for different values of T 
and V for the LM, C and MCK fixed points. The grey 
curves are obtained by varying T at fixed V for the lowest 
value of V considered in our study, i.e. U min = 10 ~ S D. 
The temperature dependence of the linear response con¬ 
ductance is reproduced at all fixed points when the non¬ 
linear conductance is taken as a function of T e g. This 
is true even for values of V several orders of magnitude 
larger than V min . 

In conclusion , we have addressed the steady-state dy¬ 
namics near unconventional quantum criticality. We 
found that in the scaling regime of all the fixed points 
considered, all observables studied (y, <j> s ,G) scale in 
terms of the same but fixed point specific effective tem¬ 
perature T efI . The local spin-spin correlation function 
X and the singlet-strength cf> s assume their equilibrium 
scaling functions even far from equilibrium when scaled 
in terms of T eff , i. e. T eff replacing T. A similar result re¬ 
lates the linear and non-linear conductance. We note that 
in the (non-interacting) pseudogap resonant level model 
such behavior is absent |46| . It has been shown that the 
non-equilibrium current noise near quantum criticality in 
models possessing gravity duals appears thermal [35, 35] , 
Our results imply that similar results hold for a larger 
class of quantum critical systems and quantities. The 
results reported here may thus help in identifying uni¬ 
versality classes of unconventional quantum criticality. 
To which extend our results rely on locality needs to be 
further investigated. 
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GENERATING FUNCTIONAL ON THE KELDYSH CONTOUR 


The generating functional on the Keldysh contour can be written as 

Z[£] = J Dcj Df J /JAe^ 9 ' 7 ’ c+/ V /) e iqN A dzA * 

Xg +i W / T dz /cL C 0,a<W'*)(XL' elc'ljv'*) 

xe cHc+atc+fU f +z}f ( S .i) 

where £ c and act as sources to the fermionic c and / fields and A is a scalar Lagrange multiplier enforcing the 
constraint Yta ft ft? = 9 -W- f 7 dz is the integral over the Keldysh contour 7 = 7_> + with its forward (7^) and 
backward (7<_) branches. 

Here, the inverse bare propagators are 

gj 1 = ( id z - \ z ), (S.2) 

gf 1 = {idz - s p i ). (S.3) 

In analogy to the equilibrium procedure -albeit performed on the Matsubara contour- one can introduce a Hubbard- 
tratc 
Thus, 

Z[£]= J Df J DX J DBe i { f ' 9 f lf )e i f dzX * Q e i { B « 9 B lB ‘f) 

Xe "4 / dzf dz' f2 aal B al zBY z ,ge,l{z,z')fl,Jvz' 


Stratonovich decoupling field B a i z conjugated to Yta' c o aa'i-f&'zt to decouple the quartic fermionic term in Eq.(S.l I. 


Xg[ yW A ’^-•x'a'{lX-9o) 0a , a , l , z B] xllz ,f„i > .i Jjj f y dzJ2^ a fl , z B a i z (g o .£ a ) 0arrlz \ 
Xe *^9cScg/ t U+?// 


with 


where 


Ob-w — 


J 


-1 


J w 


J 


= Jip, is the bare inverse propagator of the B field. 


Finally, with the help of the complex-valued dynamic Hubbard-Stratonovich fields W zz i-i one obtains 

Z[£]= J DX j DWeNtrlnl-iGJ^-Mtrlnl-^+VlGfV^)] 


iJVtr[w t *[g c ]- 1 *w]+* f dzA z Q 
e -ii)G,V ic [G^ 1 +VlGfV ic \- 1 vlG t £, f +i(,lg c i a +i4} J G f (, f 


(5.4) 

(5.5) 


(S.6) 


Gf{z,z') = g f (z,z') -W zz t ;( 
Gb 1 (z, z') = g^iz, z’) - W z t z g 


tr 


W f * [g c 


*W 


E 

1 ■ 


dz / dz 


,W zz '.iW zz 'g 

9c,i(z,z') 


(5.7) 

(5.8) 

(5.9) 


with 
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and and V^ c are source-dependent terms. Eq.(S.6) is used to derive all correlators by taking derivatives with 
respect to the source fields. 


DYNAMICAL LARGE-N SELF-CONSISTENCY EQUATIONS ON THE KELDYSH CONTOUR 

In this section we set the sources to zero and compute the saddle-point equations with respect to the bosonic fields 
W and A. The generating functional in the absence of sources is 


Z[i = 0] = J DX j DWe iNS[W ' X] 


with 


S [W, A ]=q dz\ z + tr * [g c ] 1 * W 


— i —tr In 
N 


-ZG7 1 


+ in —-trln T—zGr 1 ! . 
M 1 B J 


(S.10) 


(S.ll) 


The saddle point equations are obtained by putting the linear variation of S [W, A] with respect to W and A to zero: 

V,\\ = W zz ,, l [g c {z,z')]~ 1 (S-12) 

a 

S [W, A] = [g c ,i{ z ,z')]~ 1 W zz - ! i - iK^- ^G B ,u{ z ,z') = 0, 


SW ZZ v 

sw: 


[W, A] = q + Y, G f( z ~’ z ) = °- 


These equations become exact in the large-N limit. These equations are equivalent to 

Gf 1 = gf 1 -Y,f 

G b = 9b 1 “ ^ b 
q = —iG f (z~,z) 


with 


y /_ „/n { W Z ' Z>L 0 \_ ( g CtL (z',z) 0 

Eb(2 ’ 2) -1 o 0 

Zf(z,z') = = S aa riK y; g c> i(z, z')G b -u Q, z'). 


Gf{z,z') 


(5.13) 

(5.14) 


(5.15) 

(5.16) 


Note that A^ evaluated at the saddle-point is time independent, i.e. A t = A. 


SINGULAR EXCHANGE COUPLING MATRIX J 

So far, the treatment has been general and no particular form of the Kondo exchange coupling matrix has been 
assumed. For the physically most relevant case where the Kondo Hamiltonian is derived from an Anderson-type model 
through a Schrieffer-Wolff transformation, the exchange matrix Jw (1,1' = L,R) takes the from 

j_( Jl VJlJr\ 

V V JlJr Jr ) 

Thus, the exchange coupling matrix is singular, det( J) = 0. In this case, where one of the eigenvalues of J vanishes, 
we can write 

j = \u+)(J l +Jr)(u+\ 


with 



















As the exchange matrix is singular, the component u_ of the B field has to vanish and thus 


Gb = |w+) Gb+ (u+| . 


In this case the self-consistent equations simplify to 

E b +(z,z') = -ig c , + (z' ,z)G f (z,z'), 
Zf(z, z') = iK,g Ci+ (z, z')G b ,+ {z, z'), 

with 


G b + ~ 9 b+ 

— 1 _ ~1 

9b+ ~ Jl + Jr 

_ JL~g c ,L + Jr<]c,R 
9c,+ - Jl + Jr 

Using Langreth’s rules, we obtain 

E>’< 0 t,t l ) = -i~g<’>V,t)G>’ < (t,t l ), 
(t,t') = iK~g>’< 
q = —iGf (0). 


(5.17) 

(5.18) 


(5.19) 

(5.20) 

(5.21) 


NON-EQUILIBIRUM STEADY-STATE (NESS) DESCRIPTION 


The steady-state condition implies that the system is time translationally invariant so that G R,A,K (t, t') = 
G R,A,K (t-t'). Therefore, it is advantageous to solve the self-consistent equations in the frequency domain. The 
conventions of the Fourier transform used by us are 


A (t) = J 

A(w) = J dt A (t) e iuit . 


Eq.( S.19|S.21 1 take the form 


E>’< ^ - 


M = -i J 0 ~ w) G^’ K (u) , 

/'< M = iK J “ u ) G b < i v ) > 


(5.22) 

(5.23) 

(5.24) 


The reservoirs are in equilibrium and are thus characterized by their respective chemical potentials and pr and 
their respective temperatures Tl = Tr = T . We introduce the following reservoir quantities 


P% M = 


1 

27TJ 


Pc,l M = --P 

7r 


9c,i M ± 9c,i M 

p7,i O') 


dv- 


UJ — V 


(5.25) 

(5.26) 


where p c t (w) = ^ J2 P 8 (uj — £ p i) is the normalized (f dojp c ; (to) = 1) local density of states of reservoir l, p R t (w) 
is its Hilbert transformed and (w) is proportional to the Keldysh component of the Green’s function. Since the 
reservoirs are taken to be in equilibrium, the fluctuation dissipation theorem can be applied and it is found that 

Pt,i M = f i M Pc,i M > 


(S.27) 
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with 


f i M 


[1 — 2 iif t (lo — Hi)] = tanh 




(S.28) 


Here, n /, (w — hi) is the Fermi-function, and pi and hi are the inverse temperature and the chemical potential of 
reservoir l. The lead’s Green’s functions can thus be written in the form 


with 


9c’ A M = -7T [Pc M ± % M] . (S.29) 

9c M = -27rip+ (w), (S.30) 


Pc M 

P^ M 


j lP c ,l M + M 
Jl + Jr 

J L i L (uj) p~ L (uj) + (u;) (uj) 

Jl + J r 


(5.31) 

(5.32) 


Self-consistent equations for the steady-state 
With the definitions of the previous sections, Dyson’s equation translates to 


Pf (uj) = oj (uj) < uj — A + 7T(jjr (u>) 


2 


2 l 

+ 

7T(J f (<jj) 



-l 


Pb M = Vb M { [- ( J L + Jr) + 71-crf (i/)] 2 + [ttct b (w)] 2 } , 

with A = At — | (Jl + Jr) f du>p+ (uj) being a renormalized chemical potential, and Eq.(|S.22|s.24l translate to 


a B ( w ) _ T; 


dv 


pi (v - u) p + f (v) - pf(V- uj)p f (u) 


af (uj) = ni J dv [pf (uj - v) p% (u) + pt (uj - v) p B (v)\ , 

1 — J dujp~f (uj) . 


q= 2 


(5.33) 

(5.34) 

(5.35) 


In the particle-hole symmetric case (q = 1/2) and for a particle-hole symmetric DOS of the leads (p c (uj) = p c (—uj)) 
the quantities p± B and cr jr B are real. 


Details of the numerical treatment 


The explicit form of the pseudogap density of states of the leads is taken to be 


Pc,i M 


V2AT (r±i) 


V2A 


r 

e 


2 



with l = R,L and A = 1 specifies the soft high-energy cutoff. The self-consistent equations were solved iteratively 
on a logarithmically discretized grid with 350 points ranging from — 10A to 10A. The criterium for convergence of 
the selfconsistency loop was that the relative difference of two consecutive iterations was less than ICG 6 . The results 
were benchmarked by the conditions that the fluctuation dissipation ratios of the Green’s functions have to accurately 
reproduce the equilibrium fluctuation dissipation relations demanded by the fluctuation-dissipation theorem. For 
all the fixed points we studied a range of temperatures T/D = 10~ 1 ,10 -15 ,10 -2 ,..., 10~ 8 and a range of voltages 
T/D = 10 —2 ,10~ 15 ,10~ 2 ,..., ICG 8 . However convergence of the numerical solution of the self-consistent equations 
was not always achieved for all combinations of parameters. 
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OBSERVABLES 


Cross ^-point function 

In order to compute the currents and the Kondo singlet strength we will need to evaluate the connected 4-point 
function (T'y c p 2a2<T2 i 2 c Pl <*, a-, I, ft 2 fs^ ^ • Here, C denotes the connected part of a correlation function and T 7 is the 
time-ordering operator on the Keldysh contour. Using the procedure outlined above, one obtains 

1 1 


(T 1 f Sl (U) f\ 2 (t 2 ) 

c piaiaih (h) 

C P 2 a 2 <T 2 h Q N n ^ Sia2 ^ S2ai ^ aia2 ^P lll ’P 2lj2 


and 


^pih;P2l2 (UU 2 U 3 U 4 ) 


y/JfoJi 1 


i/( Jl + </«) (</l + Jr) 
x I dZ ' I dZ Gf 9c ’ P2h ( 2 ^4)G'b {z,z') g c - Plh ( t 3 ,z')G f (z',t 2 ) 

with 5c;pPi (i, t 7 ) = (fpUi| ffc I^PiU)- For equal times we have 

{ C P2a2CT2i2 (^) C Piai<Tlil W /s 2 (0 /s 1 c = (t) ) 

where the time-ordering for the equal-time limit is defined through F pih;p 2 i 2 (f) = 

li m *i, 2 , 3 , 4 ^t ipp^pais. (* 1 ,* 2 ,*3,* 4 )1 1 >t >t >t • Fpih;p 2 l 2 (t) can be explicitly evaluated using Langreth rules 
and making use of the fact that we describe a steady-state. This procedure is straightforward but involved and yields 


F Pih;p 2 h (t) = 4 i 7 T 5 


r(b 


hPi,l 2 P 2 I 1 P 1 J 2 P 2 


% 


( 2 ) 


with 


- 1 p^l) / % < H* IV + 1 <«) + -V + - (»)] [2p£(")l [«[ V + 1 <«) + V M] 

+ l-m [Hr 2 ++-] (w) + AJ-++- (w)] [-ipf (w) + Pb (w)] [ift [A7'++] M + H ; “ ++ (w)] 

+ [-m [4” ++ ] M + H--++ (w)] [ipf (w) + Ps (w)] [ift [H-++-] (w) + H-++- (w)] } 

J % {[-^ [V + 1 m+ai ++ “ m] at ++ m 

- [-^ [AT ++ ] M + -4, 2 V • 1 M] A-++- (w)} , 
where we defined 

A ?(u) = j 7^ P/ (1) iy)p^n (y-u) -P/ (3) {v-u) , 

U[A](u) = --P [ dv A ^ 
n J c 


OJ — V 


Currents 

The currents of particles and energy through the system are obtained from the change in particle number and 
energy of e.g. the left lead through a continuity equation for the conserved charge (particle number or energy), 

J b = ~d t ( Q b (t)) = -i ([H ( t ), Q b (t)}) 

Je -a Qe = H l = £ pr/- paaL Cp anL 

p,ot.cr 

J\> -> Qp = N L = ^ cl aa L C pacrL 

P,Oi<7 
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Using the identity [c^cg, cj^ca] = Sp^c^cs — 8 a ^cLcp and the fact that the Hamiltonian can be decomposed as 
H = Hl + Hu + Hj with 

Hj = — ^ Jw ^ ^ (/<J fa' — qSaa') c o : a<r'l' C Oaal, 

IV gc r' Oi 

one obtains 


Jv (t) /M = 2 yj Jl, 4 J_R j4 Re 


■2 Fr P', L P (*) 

'‘'C / 

pp 


Je (T) /M = 2 JL,tJn,tRe 



^ \ BpL^Rp' y p if) 
pp' 


Susceptibility 


On the Keldysh contour the impurity spin susceptibility is defined by 


(^') = -4E<V° (z)S a (z')), 


XI 


IV 


where T 7 is the time-ordering operator on the Keldysh contour. 

For a steady state, we obtain 

x -± ( w ) = J dv p+ {u - w) pf {v) - p~ f (w - u) pj (y) 
where X f M =-^ [x> H ± X< M]. 

Kondo singlet strength 

It follows from the Hamiltonian, Eq. 0> that the Kondo term contribution to the total energy is given by 


< s w - s ^' (*)> 


= — Jk 


W c 

N 2 - 1 
N 

/TV 2 - 1 

V N 


^ -hi l, 

lih 

(j) s (t). 


2 ^ / Fp 1 l 1 - I p 2 l 2 ( t) 


This expression can be greatly simplified using the definition of F Pi i i;P2 i 2 (t), see previous section. This then yields 
for the Kondo singlet strength 


with 


(j> s = n/J w h (0), 

W~ (u) = \ J dv [ a B 0 - w ) Pb 0) ~ a B^- w ) Pb M] ■ 

ADDITIONAL NUMERICAL RESULTS - OTHER VALUES OF r AND k 


In this section we provide further numerical support for our conclusions. Figure |S1| shows our results for the 
parameter set (r, k) = (0.15,0.16) which is different from the one the results in the paper are based on. 
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x'(o)- 1 /? (a) 

10 ° 

10 -' 

1(T 2 

10~ 3 10- 2 10 -' 





T e «■/D Toff/D 


Figure SI. Scaling of different observables with Tj e g for the different fixed points (the parameters used here differ from those 
of Figure 3 and 4 of the paper): (a) Inverse static susceptibility \' (0) _1 vs T e g ; (b) di nul lny" (u>) vs to/T e g; (c) singlet strength 
4> s vs Tiff For each fixed point, the equilibrium scaling form (black dashed lines) is compared with the same quantity under 
non-equilibrium conditions where T is substituted by T e ft. (d) Conductance G as a function of temperature computed for the 
lowest non-zero value of V for several values of J (see color coding), (e) G = Jv/V vs Teff- for the different fixed points. The 
equilibrium form is depicted by the black dashed lines. Go is defined as the zero-temperature limit of G in the MCK regime. 
























